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Coherent control and the creation of entangled states are discussed in a system of two super- 
conducting flux qubits interacting with each other through their mutual inductance and identically 
coupling to a reservoir of harmonic oscillators. We present different schemes using continuous-wave 
control fields or Stark-chirped rapid adiabatic passages, both of which rely on a dynamic control 
of the qubit transition frequencies via the external bias flux in order to maximize the fidelity of 
the target states. For comparison, also special area pulse schemes are discussed. The qubits are 
operated around the optimum point, and decoherence is modelled via a bath of harmonic oscillators. 
As our main result, we achieve controlled robust creation of different Bell states consisting of the 
collective ground and excited state of the two-qubit system. 
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I. INTRODUCTION 

Superconducting solid-state qubits are promising can- 
didates for quantum computation and quantum informa- 
tion, for example, because of their inherent scalabihty 
using well-established micro fabrication techniques and 
the ability to design them to meet specific characteris- 
ticsii^. The entanglement of multiple qubits is at the 
heart of quantum computation and quantum informa- 
tion, and thus has been studied extensively in the past 
few years. For example, entanglement between a super- 
conducting flux qubit and a quantum harmonic oscillator 
such as a superconducting quantum interference device 
(SQUID)^ or a LC circuil4 were examined, as well as be- 
tween two qubitsi. Entanglement or coupling between 
two superconducting charge or flux qubits are well stud- 
ied in theorji^i^ and also observed in experiments^i^i^ii^. 
Using Josephson charge qubits coupled by an induc- 
tance, a scalable quantum computing architecture was 
proposcdii. Recently, a superconducting quantum de- 
vice consisting of four coupled flux qubits was achieved 
experiment allyi^ . 

A key limiting factor for all of these devices is deco- 
herence such as depliasing or energy relaxation, which 
occurs, e.g., when the quantum devices couple to envi- 
ronmental degrees of freedom at a flnite temperature. For 
example, charge qubits are very sensitive to background 
charge fluctuationsiii^. Flux qubits are practically insen- 
sitive to background charge fluctuations, but their phase 
coherence can still be destroyed by a large number of ef- 
fectsi^. Therefore, although superconducting qubits with 
long decoherence times up to the order of jis are avail- 
able^>^'^>^>^, applications are severely limited by the 
relaxation time. Zhang et al. proposed a method to pro- 
tect stationary entanglement in superconducting qubits 
from the relaxation and dephasing processes^^. They 
found a maximum stationary concurrence of about 1/3 
and fidelity of 2/3. 

Another issue is the realization of robust quantum op- 
erations. To date, operations on SQs are typically based 



on special-area pulse9^ i^°'^^i^° , which require an accurate 
control of field parameters such as intensity, duration 
and shape. In a recent experiment, pulse-timing uncer- 
tainty limited the achieved fidelity by an error of about 
10%ii^. Similar problems are faced in the preparation 
of atomic systems, where coherent population transfer 
schemes have been developed in order to overcome these 
limitations2^i2^. In contrast to atoms, superconducting 
qubits have the advantage that the transition frequen- 
cies can be changed to a large degree on demand via the 
bias flu5*^'2i'i^, bias charge^ or bias current^!. Making use 
of this advantage, very recently, based on the breaking 
of parity symmetries, a coherent population transfer in 
superconducting current-biased phase qubit was demon- 
strated using so-called Stark-chirped rapid adiabatic pas- 
sages (SCRAPs)2^i2ii25. 

Here, we discuss schemes to create a set of relevant 
collective states in a system of two inductively coupled 
flux qubits driven by time-dependent magnetic fluxes 
(TDMFs) , addressing both the questions of robust coher- 
ent control and of decoherence. Controlled state transfer 
is reported using three different techniques: Continuous- 
wave control fields, robust SCRAP-based state transfer, 
and special-area pulses. A dynamic control of the qubit 
transition frequencies via the external bias flux is applied 
in order to maximize the fidelity of our target states. 
The inevitable decoherence is modelled via an interac- 
tion with a reservoir consisting of an ensemble of har- 
monic oscillators. Our approach has the advantage that 
the qubits operate around the optimum point, where the 
energy levels are symmetric as a function of the bias flux. 
In addition, flux qubits typically feature longer decoher- 
ence times compared to current-biased qubits^^ii^. 

In particular, we first demonstrate how continuous- 
wave driving fields can be used to populate the antisym- 
metric collective state for a pair of identical qubits, which 
is usually decoupled from external electromagnetic field 
and thus hard to populate. Second, we demonstrate effi- 
cient population of the collective symmetric state, both 
directly and via the anti-symmetric state. Finally, we dis- 
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FIG. 1: Two superconducting flux qubits interacting with 
each other through their mutual inductance M and damped 
to a common reservoir modeled as an LC circuit. The indi- 
vidual bias fluxes are varied dynamically in order to control 
the qubit transition frequencies around the optimum point. 
In the figure, crosses indicate Josephson junctions, whereas 
the bottom circuit loop visualizes the bath. 
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The gauge-invariant phase drops across the three junc- 
tions in the Ith qubit are 0^''' , (/)2 and ' . Both qubits 
experience the same TDMF 

<Peit) ^ A cosiuct) , (2) 

but an individual bias magnetic flux ^''p is applied 
through each qubit. 

The qubits can be visualized as artificial two-level 
atoms coupling to a common reservoir of quantum os- 
cillators. Using the phase constraint condition through 
the ^th qubit loop 
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+ (2^$i'V*o) + 27r$e(i)/$o = , (3) 

1=1 

where $o = h/{2e) is the flux quantum, the system in 
Fig. [T] can be described by the total Hamiltonian H = 
Hq + Hb- The two-qubit Hamiltonian Hq in two-level 
approximation and rotating wave approximation is given 
byia 



cuss controlled population of the collective excited state 
and of different Bell states composed of the collective 
ground and excited states using SCRAP. Interestingly, 
in contrast to atomic systems, the population of the col- 
lective excited state from the ground state is possible di- 
rectly without involving the intermediate symmetric and 
anti-symmetric states. 

The outline is as follows: In Sec. HIl we introduce the 
Hamiltonian of our system, which is general in the sense 
that different mutually exclusive processes are modelled, 
which become relevant depending on the parameters of 
the applied fields. This Hamiltonian will successively be 
used throughout the different control schemes presented. 
In Sec. IIIH the antisymmetric state is populated with a 
near unit population for two identical fiux qubits. Sec.lTVI 
studies the preparation of the symmetric state both from 
the ground and from the antisymmetric state. In Sec- 
tion |Vl we propose a robust method to create the col- 
lective excited state and superpositions of different Bell 
states. Finally, Section |Vl] briefly discussed and summa- 
rizes the results. 



II. MODEL 

Our model system consists of two flux qubits coupled 
to each other through their mutual inductance M— and 
to a reservoir of harmonic oscillators modeled as an LC 
circuit, see Fig. [T] Each qubit loop contains three junc- 
tions: two identical ones and one which is smaller by a 
factor of . The Josephson energies and capacitances in 
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The transition frequency luq"^ of qubit I is determined by 

^icjq'' = V^f ~^ ^1 with the tunnel coupling ti between 
two wells in the Ith. qubit and the energy difference ei of 
the wells measured with respect to the degeneracy point. 
This frequency can be expressed as ei = 

with the persistent supercurrent /^'^ and the bias flux '^^P 
in the Ith qubit loop^. The Pauli matrices of the Ith qubit 
with groimd state \gi) and excited state je;) are defined 
as 
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The phases and amplitudes of the coupling strengths 



ki^n^l^' and fil^l^ can be controlled by the applied 
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TDMFs. The always-on coupling parameters are given 



h\2 = M ( ei, 62 
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Since IA2I ^ w;, terms proportional to A2 and its complex 
conjugate can be neglected. Our system works near the 
optimal point 



/, = $(0/$o = l/2, 



(7) 



where the parameter Ai determined by the persistent su- 
percurrent is a real number— i^^. 

The interaction between flux qubits and the reser- 
voir can be described by the Jaynes-Cummings Hamil- 
tonianii^>25^>28^. 
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The bath oscillators have frequencies ujr, and {alar) = rir 
is the average photon number of the rth field mode. The 
thermal average number of photons in the oscillator 



Nth{uj) = [cxp{hw/kBT) - 



(9) 



is assumed negligible at the frequencies relevant to our 
system. Tyr'"* and Xr ™'' are determined by the vacuum 
field. In our model, we dropped the vacuum-induced de- 



cay with difference frequency \ujq 
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. Since the mag- 



nitude of the vacuum field is proportional to the square 
root of its frequency's, and because jwg^^ — Wg^^j <C 

Lo'^\ Wq^^ -t-WQ^"*}, processes at the difference frequency 
are highly suppressed. 

The diagonalization of the Hamiltonian in Eqs. ([4]) and 
([5]) without TDMF leads to the eigenenergies Ej and 
corresponding eigenstates \ j) (j € {e, s, a, g}) ^^^^^ , which 
can be interpreted as a single four-level system, 
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A = Ai — rirf includes the always-on coupling Ai and the 
shifts induced by the bath-induced dipole-dipole interac- 
tion (DDI) between the two qubits ^Id- For A = 0, the 
coefficients in Eq. ([TU]) evaluate to a = /3 = 1/ ^/2. 

Applying the Born-Markov approximation to eliminate 
the bath2^i2£i2i, we describe the bath-induced dissipation 
by the Liouville operators and ^rp, where p is the 
density matrix. The first part Ji'p results from the term 



proportional to 77^''' in Eq, 
terms2^i^ 



It is composed of three 



dp 
dt 



(12) 



where 



dp 

dt 



-£] =-Ts{{Ree+Rss)p + p{Ree+Rss) 



dp 
dt 



— 2{RsePRes + RgsPRsg)} 

+ (2a/37o + ll2){RsepRsg + RgsPRes) , (13a) 
= -Ta{{Ree + Raa) P + p{Ree + Raa) 

— 2{RaepRea + RgaPRag)} 

— (2a/37o - ll2)[RaePRag + RgapRea) , (13b) 



dp 
'dt 
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— '^{RgaPRsg + RgsPRag + RsePRea 
+ RaepRes)} + {a^ " l3^)jo{RaepRsg 
+ RgsPRea + RsePRag + RgaPRes} , 

with the damping coefficients 
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70i 7i2 are the Einstein A coefficient and the dipole-dipole 
cross damping rate, respectively. The collective qubit 
operators are defined as Rij — \i){i\ where the collective 
states |j) S {e, 5, s, a}) are given by Eq. (fTO|) . The 
contribution with subindex s [a] describes the sponta- 
neous decay via the symmetric state \s) [anti-symmetric 
state I a)]. The part with index / is an interference part 
involving both the symmetric and the anti-symmetric 
states. It results from spontaneously induced coherences 
between the symmetric and antisymmetric transitions, 
and only contributes if ^ 0^ ■ 

The second incoherent contribution J^tP takes the 
form 



■^TP = ^Qi2RgepReg - ReeP - pRee) ■ 



(15) 



It arises from the contributions proportional to Xr in 
Eq. ([5]). Note that the decay rate 70 is smaller than 
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7o according to Liu's workiS, but their ratio is tunable 
during the fabrication of a superconducting circuit. In 
order to simphfy the discussion, we assume that 70 =70- 
By applying TDMFs with different frequencies, selec- 
tive processes described by Hq become resonant, such 
that various types of effective interaction Hamiltonian 
can be generated. In the following investigation, the 

qubits work at two slowly-varying frequencies Wq^^ and 

uj'^^ , respectively, controlled by the bias fluxes The 
frequency difference A is kept much smaller than the 
always-on coupling strength Ai. As a consequence, the 
slowly-varying bias magnetic flux does not excite un- 
wanted transitions between the symmetric and antisym- 
metric states. 

Throughout our investigations below, we assume flux 
qubits with dephasing times = 1 ~ 10/is, as observed 
in recent experiments^ i^^'^^ , and an energy relaxation 
time Te of roughly half the dephasing time Tj^. Thus, 
if the bias flux is changed in order to tune the transition 
frequencies, they vary in a range smaller than 60 MHz. 
Our numerical results show that the influence of this vari- 
ation on the circulating current and subsequently 
on Ai, is negligible. Experimentally, always-on coupling 
strengths |Ai| of several hundred MHz or even higher 
are realized^i^iii^. Note that the sign of Ai can be con- 
trolled by choosing ferromagnetic or anti-ferromagnetic 
couplingi2ii2&. Therefore, it is possible to compensate the 
energy shift from the DDI via Ai. 



III. PREPARATION OF THE 
ANTISYMMETRIC STATE 

We first aim at populating the so-called antisymmet- 
ric state of the two-qubit system. For this, the frequency 
of the applied TDMF is chosen close to the average fre- 
quency ujQ. Then, the terms proportional to f^J^, 
can be dropped in a rotating wave approximation (RWA) 
in Hamiltonian 

We apply the master equation approach to the dynam- 
ics of the system in a frame rotating with the frequency 
of the TDMF ujc- In the collective states basis, the mas- 
ter equation for the system density matrix p takes the 
form 

p=l[p,Ho + Hi]+^p + J^rP. (16) 

where 

i/o ^n[5{2R,, + Rss + Rao) + w{Rss - Rao)] , (17a) 
Hi =h{a + P) [0(i?e. + Rsg) + H.c.)] 

+ h{a - f3) [n{R,a + Rag) + H.C.] . (17b) 

The detuning S = loq ~ lUc and the Rabi frequency Q = 
— fci = — fc2, see Eq. dH). We set the energy of level \g) to 
zero. The amplitude of 17 is f^o, whereas the phase 99 of 
does not influence the population, and thus is ignored in 




FIG. 2: (Color online) Time evolution of the population 
(solid red line) in the antisymmetric state \a). The idea is 
to prepare the anti-symmetric state while the two qubits are 
non-degenerate, and only afterwards render the two qubits 
degenerate. For this, the parameters are chosen such that 
712 = 0.998670, A = 50 70, 5 = 50 70, ^0 = 50 70. The two 
qubit transition frequencies are adjusted via time-dependent 
bias fluxes, such that the frequency difference A{t) (dashed 
black line) changes from I870 to zero as a cosine function 
during the time period 1207,^^ to 1607^^. The driving field 
(dash-dotted blue line) is turned off from its initial value flo 
in the period 1657^"^ to nS^g^. 



the following. The dissipative part ^p + I£yP is defined 
in Eqs. (HH) and p5)) . 

The antisymmetric state is of interest since it is par- 
tially decoupled from the interaction with the reservoir as 
Y s < 7o/2 and therefore more stable against decoherence. 
However, as is well known from atomic systems, this de- 
coupling at the same time makes a controlled population 
of this state difficult, as it also decouples the state from 
driving fields. This in particular holds if the two involved 
qubit transition frequencies are identical. 

But as an essential difference of the flux qubit system 
to atomic systems, the transition frequency of a flux qubit 
can be individually controlled by its bias flu5*^'2i'i^. In the 
following, we exploit this feature and demonstrate that 
the antisymmetric state can be populated even for identi- 
cal flux qubits by flrst adjusting the bias fluxes such that 
the qubits become unequal, then preparing the antisym- 
metric state, and finally switching back to the degenerate 
case. 

An example for this is shown in Fig. [2l Initially, the 
two qubits have a frequency difference A(t = 0) = Aq = 
I870. Applying a continuous TDMF during < 7oi < 
165 allows to populate the antisymmetric state, as can 
be seen in Fig.[2l After a certain time (7ot = 120 in our 
example), the bias fluxes are continuously adjusted such 
that the two qubits become degenerate, A(7oi > 160) = 
0. It can be seen from Fig. [2] that a preparation fidelity 
for the antisymmetric state in the degenerate two-qubit 
system of about F = 0.94 is achieved. Finally, the TDMF 
is switched off as well in the time period 165 < "fot ^ 
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175, demonstrating that it is not required to preserve 
the population in the antisymmetric state. It should be 
noted that this scheme does not rely on a delicate choice 
and control of parameters, as it is the case, e.g., for state 
preparation via special- area pulses. 

The limited preparation fidelity and the slow decay of 
the state is due to the fact that the decoupling of the anti- 
symmetric state is not perfect. For any realistic system, 
the distance between the two qubits remains finite, such 
that Ta > 0. For our parameters, 712 « 0. 998670, and 
the decay rate of the antisymmetric state is found to 
be Ta ~ 0.001470, strongly suppressed by the dipole- 
dipole coupling. In this example, the energy shift on 
the antisymmetric and symmetric states induced by the 
direct DDI is of order 10'^ 70. As discussed before, it can 
be compensated via Ai to yield a relatively small A. We 
found that the maximum fidelity obtained in our example 
is insensitive to A in a range of about A = 4O70 ~ IOO70. 

Finally, we note that the antisymmetric state is also 
an entangled state. We use the concurrence G^^^^ as an 
entanglement measure, which is given by 

C=2max{0,x(t)}, (18a) 

X{t) = \aP{pss - Paa) + a^Psa - 0^ Pas\ 

- ^/PeePgcj ■ (18b) 

In our example, C approaches 0.89 at time t — 1607q"^. 



IV. PREPARATION OF THE SYMMETRIC 
STATE 

In this section, we discuss the preparation of the sym- 
metric entangled state \s). Three approaches are com- 
pared. First, we directly prepare the symmetric state 
from the ground state using a special-area pulse. Sec- 
ond, SCRAP is used to populate the symmetric state 
from the ground state. Finally, \s) is populated via the 
anti-symmetric state discussed in the previous scction lllll 
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FIG. 3: (Color online) (a) Time dependent population of the 
symmetric state |s) directly from the ground state. The dy- 
namics is induced by a continuous TDMF with parameters 
Q.0 = 1570 and 5 = — 5O70, and the other parameters are 
712 = 0.998670, A = 50 70, and A = 0. The numerical result 
(solid red line) is well fitted by Eq. (|19|l shown as the black 
dash-dotted line. A small part of the population is trans- 
fered to the collective excited state (dashed green line), (b) 
Corresponding results for the concurrence. 



A. Direct population via special-area pulse 

First, we populate the symmetric state directly from 
the ground state. If the two qubits are closely-spaced, 
the dipole-dipole level shift is large enough for an almost 
selective excitation of the collective states without pop- 
ulating the collective excited state |e). For example, it 
was shown that in such a system of two identical two-level 
atoms, the selective population of \s) is possible through 
a short 7r-area pulse, even though there are competing 
channels such as exciting the atoms from the symmetric 
channel to the collective excited stated. In the follow- 
ing, we also consider two identical flux qubits with equal 
transition frequencies, and assume that they experience 
the same driving TDMF. The TDMF is chosen to be near 
resonant to the \g) <-> \s) transition, i.e., 5 ~ —A. 



In analogy to Ref. [35|, in our scheme described by 
Eq. (|16p . the state \s) can also be selectively prepared 
by a standard 7r-area TDMF pulse with sufhciently large 
detuning to the \s) ^ \e) transition. When the system 
initially is in its ground state, the maximum population 
of state \s) is obtained at time t = 7r/(2\/2ilo)- 



In Fig. 3(a) 



we show results for the population of \s) 
for continuous driving. It can be seen that the population 
reaches a maximum value, but afterwards exhibits rapid 
oscillations at frequency 2v^rio, while the amplitude of 
the subsequent maxima in the population decays as an 
exponential function exp[— (70 + 712)^] until the system 
approaches its stationary state. This result can be un- 
derstood by reducing the system to a two-state system 
only involving in the states \s) and \g). The numerical 
results can be well fitted by the solution of this two-state 
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approximation, 



„2-1ovg1 



I _ g-(7o+7l2)t 



+ e-^'^o+'^i^)* sm^{V2not) , (19) 



as can be seen from Fig. 3(a)[ The time evolution of the 
concurrence shown in Fig. |3(b) exhibits osciUations along 
with the population of \s). As expected, the maximum 
concurrence of C = 0.83 occurs at time t: / {2^/2rio) . 

We chose to display the result for a continuous driv- 
ing field rather than for a 7r-pulse in order to illustrate 
that the maximum population of state | s) is strongly de- 
pendent on the parameters of the driving TDMF. If the 
Rabi frequency Qq or the detuning 6 are not precisely 
controlled in an experiment, then the optimum 7r-area 
pulse is not applied, and the population of the excited 
state is strongly reduced. In our example, there is an 
optimum Rabi frequency of about 1670. For resonant 
excitation (S = —A) with this Rabi frequency, the pop- 
ulation of state |s) reaches its maximum value 0.90 at 
time 0.077(^"'^. For f2o = 570i the maximum population 
decreases to 0.84, while the maximum concurrence be- 
comes to 0.79. For Qa = 2570, the maximum population 
and concurrence are only 0.78 and 0.60, respectively. For 
a non-ideal detuning 6 — — 4O70, the population and the 
concurrence have maximum values of 0.80 and 0.64, re- 
spectively. 



B. Direct population via SCRAP 

So far, as in our previous section, in most cases special- 
area pulses have been used to create coherent superpo- 
sitions in SQa^ i^°'^^i^° . However as discussed above, this 
technique is not robust: variations in pulse area and de- 
tuning from resonance can lead to considerable loss in 
preparation fidelit}*i. In order to overcome these prob- 
lems, very recently, a so-called Stark-chirped rapid adi- 
abatic passage (SCRAP) scheme for robust population 
transfer known from atomic systems^^ has been proposed 
in superconducting current-biased qubits^i. 

The SCRAP technique builds on the well-known 
method of rapid adiabatic passage and works as follows. 
First, a pump laser pulse tuned slightly away from res- 
onance with the transition between two bound states is 
applied, and then a second delayed Stark pulse sweeps 
the bound states through the resonance by inducing a 
dynamic Stark shift. Throughout this process, the popu- 
lation is adiabatically transferred between the two states. 
Also arbitrary superpositions of two bound states can be 
achieved via the SCRAP technique^. 

SCRAP is particularly useful in our system, since a 
modification of the bias fluxes allows for a convenient 
control of the qubit transition frequencies. To apply 
a SCRAP scheme to populate the symmetric state, we 
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FIG. 4: (Color online) Robust populating the symmetric 
state from the ground state via the SCRAP technique for 
712 = 0.998670, and A = 50 70. The solid black line shows 
the population of the desired symmetric state, while the thick 
solid blue line and the dash-dotted red line are the time- 
dependent Rabi frequency and detuning required for SCRAP. 



choose Gaussian shaped pump TDMF and bias fluxes, 

n{t) = f^oe-^*-"")'/^' , (20a) 
5it) ^So~ ^oe-(*--)'/^' , (20b) 

with r!o = 3270, ^0 = 18 70, Ts = Om%\ Tp = T„ 
Ts — Tp = Ts and 5q = —60 70 corresponding to an exact 
detuning —10 70. The delayed dynamic bias flux follows 
the pump pulse. For this, the position of the symmet- 
ric collective state is adjusted in time via the bias flux. 
Because the competing channel \s) <-> |e) limits the in- 
tensity and duration of the pump pulse, the population 
(black line in Fig. 2]) in the symmetric state is limited 
to about 0.81 directly after the SCRAP preparation at 
time 1.55 7(7^. "The corresponding concurrence is 0.64. 
Even though this fidelity is rather low compared to the 
other preparation schemes presented here, the SCRAP 
approach may be an alternative if suitable and reliable 
special-area pulses to populate \s) are hard to achieve 
experimentally. 

C. Population via the anti-symmetric state 



In Sec. mil we have shown that a nearly complete pop- 
ulating of state I a) can be obtained for two identical flux 
qubits, without the need for precisely controlled TDMF 
pulses. In this subsection, we extend this scheme by 
preparing the symmetric state from the antisymmetric 
state in two identical flux qubits. In comparison with 
the direct preparation of \s) from the ground state dis- 
cussed in Sec. IIV Al the population is first transferred to 
I a) and afterwards transferred to \s). If the driving field 
is applied continuously, then the system will turn out 
to mainly oscillate between the two maximally entangled 
states I a) and |s). In contrast to simple special-area pulse 
schemes, the two Rabi frequencies driving the qubits are 
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applied with a fixed relative phase, which is possible since 
we work around the optimum point, as discussed below. 

To selectively excite the transition between state \a) 
and |s) imder the condition of A 3> |A|, the frequency ujc 
of the TDMF should be close to the level shift 2A rather 
than to the frequency difference A. In this case, a rotat- 
ing wave approximation eliminates terms proportional to 
ki and Q^f^ in the Hamiltonian Eq. ([4]), such that the con- 
tribution proportional to r2|^ determines the dynamics. 
The master equation then reads 



dp 

at 



-i[nRas + n*Rsa,p]+.^P + ^rP, (21) 



where the detuning S = w — ujc/'2. The Rabi fre- 
quency isn^ a^C - 0^Q* ^ {a^ - fi^WQ] + i3?[C] with 
C, = fl'{2 + ^21''* ■ Here, 5R[ ] and 3[ ] denote the real part 
and the imaginary part, respectively. Since — w 
for small A, a large imaginary part of is required to ef- 
ficiently drive the system. Thus, the phase of the driving 
TDMF should be close to tt/2, i.e., O = iOq. This can be 
achieved by setting the phases 6[^2 ^-^d 6*21^ of the Rabi 
frequencies Q^^l] and fl'^-^ to 7r/2 and 37r/2, rcspcctivelyi^. 
We have numerically solved the Schrodinger equation to 
evaluate fi^^ and ^^21^- We found that the two coupling 
coefficients have the same magnitude but opposite signs 
if the two bias fluxes are at symmetric positions with 
respect to the optimal point /; = 0.5. Therefore, it is 
possible to individually control the phases of the Rabi 
frequencies via the bias fluxes as well as the phases of 
the driving TDMF. 

In this conflguration, the driving TDMF is well off- 
resonant from the transitions of \g) \s) and |e) <-> |s), 
such that no such transitions are induced by the TDMF. 
Still, the population is spontaneously damped from the 
symmetric state |s) to the ground state. The equation of 
motion for the population in the ground state is obtained 
as 



dpg 



dt 



70 [1 - Pgg(O)] + 2a/37i2 [pss - Paa)- (22) 



This can be solved to give 
Paait) = 1 - [1 -Pggi 



-70* 



(23) 



+ / 2aP-fi2[Pssit')- Paa{t')]dt' 

Jo 

~l-[l-p,g(0)]e-^°*. 

Here we drop the integral term which describes a small 
and rapidly oscillating perturbation. Similarly, the time- 
dependent population in the state |s) is 

p.sit) = e-^ot ^S^ + nl^j . (24) 

In the above equations, Pgg(0) and Paa(O) denote the ini- 
tial populations in the states \g) and |a), respectively. 
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FIG. 5: (Color online) Population of the symmetric state 
|s) from the anti-symmetric state \a). The two panels show 
the time-evolution of (a) population and (b) concurrence for 
712 = 0.998670, A = 5O70, A = 0, 5 = and Oo = I570. The 
populations in the state \g) (Black dashed line) and |s) (solid 
red line) are well fitted by pgg{t) and Psa(t) (fits shown by 
green dotted lines), respectively. 



To study the fldelity of the population transfer between 
I a) and |s), in Fig. 5(a)[ we choose the anti-symmetric 
state as the initial condition. Applying a continuous- 
wave TDMF with Rabi frequency Qq = 15 70 and detun- 
ing (5 = 0, the symmetric state \s) reaches its maximum 
population of 0.90 at time 0.17,^^. After this maximum, 
the population continues to oscillate between \a) and \s) 
due to the applied field. This oscillation is damped by an 
overall decay as we include damping with a ra te 70- The 
corresponding concurrence is shown in in Fig. 5(b) The 
concurrence oscillates at twice the frequency of the pop- 
ulation oscillation, since both \s) and \a) are maximally 
entangled. The local maximum values of the entangle- 

\s) 



mcnt occur at times n7r/(2y^P + fig) where either \s) or 
I a) is occupied. 

One can improve the maximum population transfered 
to \s) by increasing the Rabi frequency, because then the 
transfer to the symmetric state is more rapid and thus 
leads to less damping throughout the transfer. This en- 
hancement is limited by the fact that the Rabi frequency 
must not become strong enough to also induce transi- 
tions involving the collective excited or ground states or 
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FIG. 6: (Color online) Direct population transfer from the 
collective ground to the collective excited state (solid red line) 
using the SCRAP technique. A Gaussian pump pulse (dashed 
black line) is followed by a delayed Gaussian time-dependent 
detuning (blue thick line). 



to break the rotating wave approximation. For example, 
a population of the symmetrie state of 0.97 is obtained 
if we use a TDMF with f^o = 50 70 and 6 = with dura- 
tion 7r/(2rio), i.e., a 7r-pulse. In this case, a concurrence 
of 0.97 is achieved. 

We also found that the population transfer is rather 
insensitive to the detuning. For example, with Qq = 
50 70 and 5 = 7"fo, the population of the symmetric state 
changes only slightly. Similarly, the scheme works well if 
only parts of the population are initially in the state \a), 
as can be seen from Eq. p4l) . 



V. PREPARATION OF BELL STATES AND 
THE COLLECTIVE EXCITED STATE 

Interacting atomic systems do not allow for an efficient 
transfer of population to the collective excited state, as 
the transfer always has to proceed via either the collective 
symmetric or the antisymmetric intermediate state. In 
the following, we show that a controlled change of the flux 
qubit properties allows to robustly generate the collective 
excited state, as well as the Bell states 



l$±> = -^(|.9i,<?2>±|ei,e2)). 



(25) 



via a mechanism not present in atomic systems. For this, 
we now assume the frequency-matching condition lUc = 

-^4^^ in Eq. Q. 
Under these conditions, a process described by the 
term proportional to Q^f^^ in Eq. ([4]) is possible that does 
not occur in atomic systems, where both qubits evolve to 
their excited states after absorbing one common photon. 
The equation of motion for the density matrix becomes 







— ^ 
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FIG. 7: (Color online) Creation of Bell states |$±) using the 
SCRAP technique from the ground state. The solid red line 
denotes population p+ in |$+), while the dashed green line 
shows population p_ in |$-). The concurrence (dash-dotted 
blue line) has a maximum value of 0.94. The black dash 
double dotted line indicates the applied SCRAP pulse Rabi 
frequency, while the thick blue line shows the time-dependent 
Stark detuning. 



with 



- h[nRe 



Rss + Raa] + hw[Rs 
+ ^*Rae] , 



Ra 



(27a) 



(26) 



where the detuning 5 = ujq — uJc/'^ and Q is the Rabi fre- 
quency. Interestingly, we find that this system consisting 
of two coupled flux qubits behaves similarly to a two- 
level system composed of the collective ground state \g) 
and the collective excited state |e). Our numerical results 
show that this system is modifications of the parameters 
A, A and 712. Therefore, in the following, we use fixed 
parameters A = 0, A = 50 70 and 712 = 0.9986 70 for 
simplicity. 

In order to achieve robust control, we use again the 
SCRAP technique to efficiently transfer population to the 
collective excited state |e) and to create superpositions 
of the collective ground and excited states. In addition 
to previous workSl, we study two coupled qubits, and 
include decoherence. 

In Fig. [SI a pump TDMF (black line) transfers the pop- 
ulation to the collective excited state during the first res- 
onant crossing. The dynamic detuning (blue thick line) 
is achieved by slowly tuning the bias flux. To suppress 
the effect of the relaxation during the preparation, the 
pump laser pulse and the dynamic variation of the bias 
flux pulse ideally should be much shorter than the col- 
lective relaxation time (470)"^. We also choose Gaussian 
shaped pump TDMF and bias fluxes given in Eq. (PO]) 
but apply different parameters: ^Iq = 180 70, Sq = 16 70, 
Sq = 3O70, Ts = 0.005 7o'\ Tp = Ts, - = Ts- 
In our example, directly after the chirping of the detun- 
ing at time 1.5l7(j'^, the population of state |e) is 0.94. 
Similarly, a so-called half-SCRAP with a small static de- 
tuning (5o ~ can be used to robustly create an arbitrary 
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FIG. 8: (Color online) Creation of superpositions of Bell 
states controlled by the static detuning So. The populations 
p± in states |<1'+) (dashed red line) and 1$-) (dash-dotted 
green line) exhibit periodic oscillations as a function of So. 
The maximum concurrence C is larger than 0.95 (solid blue 
line). 



tially from an initial maximum value of 0.94. Thus we 
find that Bell states |<&-) can be created with high fidelity 
by means of the SCRAP technique. 

Finally, we show that controlled superpositions of the 
two Bell states |$±) can be prepared using our scheme. 
For this, we focus on the state of the system directly after 
the SCRAP preparation at time l.bl^Q^. The concur- 
rence and the fidelity of the Bell states \^±) as a func- 
tion of the static detuning Sq are shown in Fig. [H It 
can be seen that one can generate arbitrary superposi- 
tions of these two Bell states |$±) with corresponding 
probabilities varying between ~ 0.03 and ~ 0.96 in com- 
bination with a near complete entanglement. Because 
a very small part of the population is in the symmet- 
ric and anti-symmetric states, the sum population in two 
Bell states only 0.99, slightly smaller than unity. The 
fidelity of |$±) sinusoidally oscillates with Sq over a pe- 
riod of about 2OO70. It should be noted that the high 
maximum fidelity of 0.96 is achieved although damping 
processes are included. 



superposition of states |e) and |g) in principle^^. 

The SCRAP technique can also be used to generate 
Bell states |$±) consisting of the collective ground and 
excited states. An example is shown in Fig. [T] For the 
chosen parameters Hq — 80 70, Sq = —870, 5*0 = 60 70, 
Ts = 0.005 7o-\ Tp ^ Ts,Ts - Tp = Ts, the fidelity of 
the Bell state |$_) is 0.94 at time 1.51%^. Since |$+) 
and are not eigenstates of the coupled systen*^, os- 
cillations occur between these two states. At the same 
time, the oscillation amplitude is damped due to the de- 
cay from |e) to The concurrence is calculated in this 
case as^ 

Cit) = 2 max{0, \peg\ - ^p^^} , (28a) 

P22 = P^Pss + aPipas + Psa) + Q."^ Paa , (28b) 
P33 = a^Pss ~ af3{pas + Psa) + 0^ Paa , (28c) 

where P22 and pas are the populations of the product 
states \e1g2) and 15162). The concurrence decreases to- 
gether with the populations of the Bell states exponen- 



VI. CONCLUSION 

In summary, robust schemes to create entangled states 
or controlled superpositions of entangled states in a sys- 
tem of two coupled flux qubits have been discussed. De- 
coherence of the qubits is included via the interaction 
of the system with a reservoir of harmonic oscillators. 
First, we have shown how the antisymmetric collective 
state can be prepared in two identical flux qubits in the 
steady state using continuous- wave driving flelds. Since 
the antisymmetric state decouples from other states, it is 
useful, e.g., for the study of Bell's inequality violationsi. 
We have also compared two different channels, direct and 
via the anti-symmetric state, to populate the symmetric 
collective state. Finally, a robust and flexible SCRAP 
technique has been demonstrated to efficiently populate 
the collective excited state. This technique enables one 
to create arbitrary superpositions of Bell states |$±) con- 
sisting of the collective ground and excited state simply 
by controlling the detuning. 
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